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ABSTRACT 

Using our Binary Population And Spectral Synthesis (BPASS) code we explore the 
effects on star-formation rate indicators of stochastically sampling the stellar initial 
mass function, adding a cluster mass dependent stellar upper-mass limit and including 
binary stars. We create synthetic spectra of young clusters and star-forming galaxies 
and compare these to observations of Ha emission from isolated clusters and the rela- 
tion between Ha and FUV emission from nearby galaxies. We find that observations 
of clusters tend to favour a purely stochastic sampling of the initial mass function 
for clusters less than lOOM©, rather than the maximum stellar mass being dependant 
on the total cluster mass. It is more difficult to determine whether the same is true 
for more massive clusters. We also find that binary stars blur some of the observa- 
tional differences that occur when a cluster-mass dependent stellar upper-mass limit 
is imposed when filling the IMF. The effect is greatest when modelling the observed 
Ha and FUV star-formation rate ratios in galaxies. This is because mass transfer and 
merging of stars owing to binary evolution creates more massive stars and stars that 
have greater mass than the initial maximum imposed on the stellar population. 
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1 INTRODUCTION 

A key problem when modelling stellar populations is how 
to determine the distribution of initial stellar masses in the 
population. The conventional method is to define an initial 
mass function (IMF) according to which the number of stars 
of a given mass is calculated as a function of the mass. Typi- 
cal ly a powe r - law o f the mass is used. The first was suggested 
by ISalpeterl j 19551 ) where dN(M) oc M~ 2m dM for 0.3 < 
M/Mq < 10. Despite being 57 years old this IMF is still 
widely used and appears to be universal. This slope holds 
over a wide range of stellar masses, only fl attening in gra- 
dient below stellar masses of around IMq llMiller fc Scald 
19791 : lKroupall200ll : IChabrier|[2003l : iBastian Covey fe Meveij 



2010) 



The IMF provides the distribution of stellar initial 
masses in a stellar population such as that found in stellar 
clusters. However when trying to simulate the stellar popu- 
lation in a galaxy it is important to recognise that a galaxy 
is not made up of one unique stellar population. A galaxy 
is actually made up of numbers of stellar clusters each with 
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their own mass and age. The masses of these clusters are also 
described by their own cluster initial mass function. Fur- 
thermore fluctuations in the IMF of these clusters, especially 
if their mass is less than 1O 4 M0, can pr ovide large variation 
in the ionising fluxes from the cluster ([Cervino et a 3 120031 : 
IVillaverde. Cervino fc Lu ridiana 2OIOaj0). The implication 
of stars forming in clusters is that the stellar IMF (SIMF) 
does not apply across an entire galaxy. Instead to get the 
galaxy-wide distribution of stellar masses we must model a 
number of stellar clusters with different masses according 
to a cluster initial mass function (CIMF) and within each 
clusters apply a SIMF to produce an integrated galaxial 
initia l m ass function (IGIMF) (e.g. Weidner fc Kroupal 



2007) 



120061 ; |H amm-Altenburg. Weidner fc Kroupal 
This is the case even if not all the clusters are 
dense eno ugh to remain bound over their life- 



times JPortegies Zwart. McMillan fc Gielesl l2010l : 
Bressert et all l2010l: iGieles fc Portegies Zwartl [201lh 



Bastian. Covey fc Meveil (|2010h and lHaas fc Anders! (|2010h 
reviewed and investigated the importance of combining a 
CIMF and a SIMF to make a galaxy- wide SIMF. They 
found that the resultant galaxy-wide SIMF is most sensitive 
to the minimum mass for a cluster, the slope of the CIMF 
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and whether the mass of the most massive star is limited 
by the cluster mass. 

The most extreme articulation of the last factor is 
whether it is physically possible for a IOOA/q cluster to be 
composed of a single lOOM© star or the more probable sce- 
nario of a cluster composed of many lower- mass stars.. The 
alternative is that such low-mass clusters can only put say 10 
percent or less of their total mass into the most massive star. 
If the former is possible it would suggest that some O stars 
might form in isolation (here taken to be in isolation with 
no other O or B stars in the same cluster so the other clus- 
ter members are of a much lower mass) and this would have 
important implications for the star-formation process. I.e. 
is star formation a pure-stochastic process or a bottom-up 
process with low-mass stars formed first and high-mass stars 
only formed if there is enough material left. After account- 
ing for runaway O stars contaminating the app arent number 
of isol ated (with no companions at all) O stars I de Wit et all 
(2005) suggested that at most 4 ± 2 per cent of O sta r s form 
outside a cluster environment. iParker fc Goodwin! (|2007T ) 
considered that an isolated O star only meant there was 
no other OB star in the same cluster. This is the case when 
a lOOAf© cluster is composed of one star that contains most 
of the mass of the cluster with a few very low-mass compan- 
ions. They modelled the populations of star clusters using 
a standard CIMF with a slope of —2 and predicted that 5 
per cent of O stars formed in clusters that had no other O 
or B stars and would be observed as isolated when in fact 
they are just massive stars that have been able to form in a 
low-mass cluster. 

The argument against the conclusion that O stars 
can form in l ow-ma ss clusters has bee n pre sented 
by IVanbeverenl Jl982h. I Weidner fc Kroupal (|2006l ) and 
IWeidner. Kroupa fc Bonnelll (|2010l ) who have suggested that 
observations indicate that, for a specific cluster mass, there 
is a maximum possible stellar mass well below the to- 
tal cluster mass. If this were the case then, from their 
model, a IOOMq cluster could not form stars with masses 
above IQMq. Therefore when a synthetic galaxy is cre- 
ated from the convolution of cluster and stellar IMFs 
there would be a dearth of the most massive stars com- 
pared with when t here are no restrictions on the maxi- 
mum stellar mass. PA amm-Altenburg. Weidner fc Kroupal 
|2007l . l2009h have suggested that there are differences in 
star-formation rate indicators when stellar populations are 
modelled by the two different IMF filling methods. How- 
ever the observations used consider the differences that oc- 
cur at low star-formation rates and thus uncertainties and 
low-number fluctuations between observed systems make it 
difficult to determine whether the maximum stellar mass 
depends o n the cluster ma s s. However the complemen tary 
studies bvlElmegreenl (120061) , [Parker fc Goodwinf (|2007l ) and 
iMaschberger fc Clarke! |200£ ) show that similar observa- 
tions indicate that there is no evidence for restrictions on 
the maximum stellar mass in clusters. 

In th is pap er we exami ne rece nt observations 
jLee et all 120091 ; lLamb et al l l20ld ; ICalzetti et all 
|2010| ) that may provide a firmer constraint on 
how nature select s the masses of stars in clusters 
and galaxi es dPflamm-Altenburg. Weidner fc Kroupal 
20071. 120091; IVillaverde. Ceryino fc Luridianal l2010bl ; 



The novel feature of this work is that we are able to demon- 
strate how binary stars alter ou r population synthesis 
predi c tions. Recent observation s llPins onncaul t fc Stanekl 
120061 ; iKobulnickv fc Frveij 120071 ; iKiminki et all l2009h in- 
dicate that the binary fraction in young massive stellar 
populations is close to one. It is therefore vital to include 
binary stars especially those that interact. In our popu- 
lations approximately two thirds of binaries interact. We 
first outline our stellar evolution models and the method 
of our spectral synthesis. We then describe our two ways 
to determine the distribution of initial masses in synthetic 
clusters and galaxies. Next we discuss the observational 
implications of varying the IMF-filling method on the Ha 
and FUV star-formation rate indicators. Finally we present 
our conclusions. 



2 NUMERICAL METHOD 

2.1 Binary population and spectral synthesis 

We have developed a novel and unique code to pro- 
duce synthetic stellar populations that include binary star s 
ijEldridge, Izzard fc Toutll2008l ; lEldridge fc Stanwavll2009l l. 
While similar codes exist our Binary Population and Spec- 
tral Synthesis (BPASS) code has three important features 
each of which set it apart from other codes and enable it 
to study stochastic effects on the IMF. First, and most im- 
portant, is the inclusion of binary evolution when modelling 
the stellar populations. The general effect of binaries is to 
cause a population of stars to look bluer at older ages than 
predicted by single-star models. Secondly, a large number of 
detailed stellar evolution models are used to create the syn- 
thetic populations rather than an approximate rapid popu- 
lation synthesis method. Thirdly, we use as many theoretical 
inputs in our synthesis with as few empirical inputs as possi- 
ble to create a completely synthetic model to compare with 
observations. 

BPASS uses approximately 15,000 detailed stellar mod- 
els calculBted_by_yT«_CainbrW code as described 
bv lEldridge. Izzard fc Tout! (120081 ). These include single star 
and binary models with initial masses between 0.5 and 
12OM and 5 and 12OM respectively. We take 120M Q 
to be our most massive star possible because of our lim- 
ited grid of binary evolution models. Above this mass the 
mass-loss rates at solar m etallicity on the main-sequence 
are high (|Vink et al.ll201ll ) and the evolutionary timescales 
of the stars vary little as the initial mass is increased fur- 
ther. We note that, owing to stars merging our binary pop- 
ulations include some single stars that have effective ini- 
tial masses of 2OOM0 and above. The minimum binary pri- 
mary mass of 5M© is selected because initially our binary 
models were specifically created to study the progenitors of 
core-collapse supernovae. The main-sequence lifetime of a 
5Mq star is lOOMyrs, which is the period we use for the 
duration of the star-burst in our constant star-formation 
models so there is no effect from low-mass binaries in 
these models. Furthermore observatio ns indicate the binary 
fract i on decreases at the l ow masses ( Duauennov fc Mayor! 



Il99ll ; iLeinert et all Il997l ; IBouv et al 



2003). 



Fumagalli, da Silva fc Krumholj l201ll ; IWeisz et all 1 201 2n 



However the 

binary fraction is more complicated than we assume here 
and is determined by a star cluster's dynamics, environ- 
ment and age. It is thought that stars of all masses can 
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form in binary or multiple systems but these can be bro- 
ken up by dynamical interactions in young clusters (e.; 
iGoodwin fe Kroupal 120051; iFregeau. Ivanova fe Rasid 12009 



We note that lHan. Pod siadlows ki fe Lvnas-Grav 



found low-mass binary stars can explain the excess UV flux 
observed in Elliptical galaxies. Such systems do not con- 
tribute strongly until lGyr after formation at which time our 
estimated UV fluxes should increase only slightly. Creating 
a new grid of low-mass detailed binary models for inclusion 
in BPASS is unnecessary for this work to demonstrate the 
importance of binary stars. 

Here we use models at solar metallicity with a metal- 
licity mass fraction of Z = 0.02. We include con- 
vective overshooting and a m ass-loss prescription that 
combi ne s the mass-loss rates of Vink. de Koter fc Lamersl 
J200l[). de Jager. Nieuwenhuiizen fc van der Hucht ( 19881 ) 
and " iNugis fc Lamersl (|2000). The binary evolution ac- 
counts for Roche-lobe overflow, common-envelope evo- 
lution, mass transfer and neutron-star kicks which 
affect the survival of binary stars after a super- 
nova. These models a re combined with the st el lar at - 
mosphere spectra of Smith. Norris fc Crowtherl ll2002h. 



lHamann. Grafener fc Liermannl ~ 20061 ) and IWestera et al.l 
[20021 ) to predict the spectra of the stellar populations. 

A significant change we make here, compared with our 
previous work, is to break from our previous assumption that 
the SIMF can be described by a simple Salpeter law over 
the entire mass range of stars. Our method requires us to 
consider that all stars are born in clusters. The mass of these 
clusters is described by a CIMF and the mass distribution 
of stars within each cluster is described by the SIMF. This 
is achieved by first picking a cluster mass and then filling 
the cluster with stars from the SIMF. We model multiple 
clusters together to create synthetic galaxies with different 
star-formation histories but with the same mean constant 
star-formation rate over a long period of time. 

We use two methods of populating the SIMF for our 
synthetic clusters. They differ by whether we limit the max- 
imum stellar mass or not. Our first method is to assume that 
any star can occur in any cluster such that, M max ^ M c \. 
I.e. the star cannot be more massive than the cluster it in- 
habits. This we refer to as pure stochastic sampling (PSS) 
of the SIMF. In this SIMF we assume a Salpeter slope of 
-2.35 between 0.5 and 120 Mq and a slope of -1.3 between 
0.1 and O.5M0. It is similar to the constrai ned sampling 
method outlined bvlWeidner fc Kroupal (|2006l ) and used by 
IVillaverde. Cervifio fc Luridiana ( 2010bl ). 

Our second case has the maximum mass of 
a star in a cluster dependent on the total mass 
of the cluster. We use the relat io n ca lculated by 
Pfl amm-Altenburg. Weidner fc Kroupal (|2007l ) which is 
given by, 



log lo (M max /M ) = 2.561og lo (M cl /M ) 
x(3.82 9 - 17 + (log lo (M cl /M )) 9 - 17 ) 



0.38, 



where M max is the maximum stellar mass possible in a clus- 
ter of mass Ma.. We therefore use M max from this equa- 
tion as the maximum mass in our initial mass function up 
to a limit of 120A/ Q in our synthetic clusters. We refer to 
this method as the cluster mass dependent maximum stellar 
mass (CMDMSM) method. The resulting clusters are sim- 



ilar to those from the sor ted-sampling method outlined by 
IWeidner fc Kroupal (|2006h . 

We note that our synthetic populations have some lim- 
itations. In Figures \T\ and [2] there are diagonal and horizon- 
tal linear features in the distribution of model populations. 
These arise at low cluster-masses and star-formation rates 
owing to the limited resolution of the stellar model initial 
masses and time bins used in our synthesis. This becomes 
most noticeable when there is only one massive star in the 
stellar population. One solution to this would be to interpo- 
late between stellar models but given that stellar evolution 
is non-linear and binary evolution is even less predictable 
we avoid spurious results from interpolations and select the 
closest model available. 

Also our binary population models are not complete and 
here we are only demonstrating the importance of including 
binary stars. For example, we do not include binaries with 
initial primary masses below 5Mq and as yet we do not con- 
sider the emission from X-ray binaries. This would provide 
another source of ionising flux that would also effect the Ha 
and UV flux ratio. The effect would be more important at 
low cluster masses and low star-formation rates where one 
X-ray binary would d ominate the entire i onising flux from 
the stellar population (|Mirabel et al.ll201ll ). 



2.2 Creating synthetic clusters 

We use the PSS and CMDMSM methods to create synthetic 
stellar populations in two regimes. In the first we consider 
individual stellar clusters with all the stars coeval. We create 
models of stellar clusters with both PSS and CMDMSM 
and investigate how they affect the Ha line flux per Mq in 
the cluster. Our process for cre ating a synthet i c clus ter to 
compare to the observations of ICalzetti et al.l ((2010) is as 
follows. 

(i) We randomly generate a cluster mass between 10 
and 10 6 Mf) from the CIMF which h as a slope of -2 
l|de Griis et al.ll2003l ; ITada fc Ladall2003l ). 

(ii) We fill the cluster with stars, the masses of which are 
picked at random from the SIMF with the maximum stellar 
mass given by PSS or CMDMSM. 

(iii) We add stars to the cluster until the total mass 
is greater than our target cluster mass. We then consider 
whether the final cluster mass is closer to the target cluster 
mass with or without the last star added to the cluster. If 
the mass is closer without the last star we remove the last 
st ar from the cluster. This is similar to the sorted sampling 
of I Weidner fc Kroupal (|2006h and makes it less likely that 
a star can be added that is more ma ssive than the tar get 
cluster mass as in the soft sampling of lElmegreenl ([2006) . 

(iv) We randomly generate the cluster age between 1 
and 8Myr. This is to match the observed age range of 
ICalzetti et ail i|2010l ). 

(v) We calculate the Ha flux for the resultant stellar pop- 
ulation. This is done with theoretical stellar atmospheres 
and stellar mod els to predict the res u ltant total spectrum 
as described by lEldridge fc Stanwavl [2009). We calculate 
the number of ionising photons from wavelengths shortward 
of 912A and convert this to the flux of Ha by assuming 
10 11 ' 87 ionising photons give rise to lergs -1 of Ha flux. 



© 2005 RAS, MNRAS OOO.rTlfTTl 



4 John J. Eldridge 



This process is repeated for many different cluster 
masses so that we can build up a picture of how Ha 
flux varies with cluster mass f or clusters aged between 1 
and 8 Myr. ICalzetti et ail (|2010h performed an observational 
study of such clusters and provide the observed mean Ha 
flux per Mq for two different masses of clusters. We com- 
pare our models to these observed populations in Section 
3.11 In the binary population case we include a companion 
for every star that has an initial mass greater than 5Mq. We 
assign binary parameters at random from a flat initial mass 
ratio and flat distribution of the logarithm of the initial sep- 
aration using the model clo sest to the parameters from ou r 
grid of models calculated bv lEldridge. Izzard fc Tout! |2008). 
We include the mass of the companion in the total cluster 
mass. 



2.3 Synthetic galaxies 

Our second set of population models are for synthetic galax- 
ies with an assumed constant star-formation rate. Rather 
than fill up the population of a galaxy according to a galaxy- 
wide IMF we create the galaxy from a set of clusters that 
each have their own individual age and stellar population. To 
create a galaxy we first pick a star-formation rate between 
10 -5 to lOMQyr -1 . We then create the synthetic galaxy as 
follows. 

(i) We pick a cluster mass at random from a CIMF which 
has a slope of —2 between 50 and 1O 6 M0. 

(ii) We fill the cluster with a stellar population as de- 
scribed in Section 2.2 and aged to between and 100 Myr, 
chosen at random from a uniform distribution. 

(iii) We continue this process until the total mass created 
in the galaxy over 100 Myr gives the required star-formation 
rate. 

(iv) With this stellar population we calculate the number 
of ionising photons from wavelengths shortward of 912 A and 
convert this to the flux of Ha by assuming 10 11 ' 87 ionising 
photons give rise to lergs -1 of Ha flux. We also calculate 
the UV flux density at a wavelength of 1500A. 

(v) From these Ha and UV fluxes we calculate an appar- 
ent star-formation rate from both and find their ratio. We 
assume a star-formation rate of lA^yr' 1 produces a Ha 
flux of log 10 (F 1 (Ha)/ergss _1 ) = 41.1 and a UV fl ux density 
of log ln (F(1500A)/ergss~ 1 Hz' 1 ) = 27.85 as in [Kennicuttl 
(119981) . 

We perform these simulations for single and binary pop- 
ulations and for the PSS and CMDMSM methods of filling 
the IMF so that the differences can be compared. Here we 
use a different range of clus ter masses based on the sugges- 
tion of lLada fc Ladal (|2003l ) that there is a turn-over in the 
mass function of molecular clouds at around 50 Mq. We also 
only consider a period of 100 Myr because th is is of the order 
of a typical star-formation burst duration (McQ uinn et al.l 
2009). We also find that increasing the age beyond 100 Myr 
has little effect on our results because it is the typical lifetime 
of stars that contribute to the FUV. We note that, when 
used to create a synthetic galax y, our CMDMSM method i s 
based on the IGIMF method of IWeidner fc Kroupal (2006). 
However we do not limit the maximum cluster mass in a 
synthetic galaxy by the total star-formation rate as they 
do in their IGIMF method. Recent investigations of the 



CIM F sugge st that there is no such dependence (|Gielesl 
2009; Larson 2009). In this work we wish to concentrate on 
whether the maximum stellar mass depends on the cluster 
mass. We have calculated IGIMF models to see the effect 
of including such a limit and find our models are in agree- 
ment with those of IPflamm-Altenburg. Weidner fc Kroupal 
l|2007h an d IPflamm-Altenburg. Weidner fc Kroupal (1 2009^ 
Also like iFumagalli. da Silva fc Krumholz ( 201lf ) we find 
that IGIMF synthetic galaxies cannot reproduce the ob- 
served spread of star-formation rate ratios. This is because 
restricting the maximum cluster mass decreases the number 
of massive stars even more dramatically than they are in our 
CMDMSM models. 

In Section 13.21 we comp are the sy n thetic populations 
to the observed galaxies of iLee et al.l (|2009l ). The novel 
feature of our approach is in not forcing clusters to form 
at the same time but allowing each to have a different 
age. This leads to much more scatter in the predicted ob- 
ser vables of our synthetic galaxies . This was also found 
by IFumagalli. da Silva fc Krumholj (|201ll ) and lWeisz et alj 
l|2012h . We have also varied the age range used for the syn- 
thetic galaxies and find that increasing the age has little 
effect on our results. Using a younger upper age limit in- 
creases the amount of Ha flux relative to the UV flux. This 
is because the stars that cause Ha emission are more mas- 
sive and typically have lifetimes of 10 Myr or less, while stars 
that contribute to the FUV continuum span a much greater 
lifetime range of up to 100 Myr. 



3 RESULTS 



3.1 The Ha from individual clusters 



ICalzetti et aD (|2010l ) suggested a novel test for determining 
how the IMF defines the population of stellar clusters. They 
studied the production of ionising photons by young clusters 
in NGC5194. If an IMF is populated purely stochastically 
then one 1O 5 M0 cluster should have the same stellar con- 
tent as a hundred lOOOM© clusters. Therefore both samples 
would have the same Ha flux per Mq of stars. However if 
the IMF of a IOOOMq cluster is devoid of massive stars due 
to a link between M max and Mci then the hundred lOOOM© 
clusters would have less Ha flux per Mq than a 1O 5 M0 
cluster. 

For our population of synthetic clusters we have calcu- 
la ted the mean Ha flux per Mq as for the observed clusters 
of ICalzetti et al l (|2010l ). Figure Q] shows our synthetic clus- 
ters as points along with the mean Ha flux. We see that at 
cluster masses below 1O 4 M0 the results diverge. With PSS 
it is possible to have one massive star making up most of the 
mass of a cluster while with CMDMSM this is not possible. 
Therefore for PSS and CMDMSM the observed mean Ha 
flux drops from the mean value of around 10 34 1 ergs _1 Mg 1 
at around 10 2 or 1O 4 M0 respectively. Therefore by measur- 
ing the Ha flux for clusters in between these key masses we 
sho uld be able to determ ine how nature fills the IMF. 

ICalzetti et ail (|2010h provide two observed values for 
their two mass bins. The first and higher value does not 
include clusters that are undetected in Ha. The second in- 
cludes these non-detections. The observations at a cluster 
mass of 10 4 ' 5 M Q agree with the predicted mean Ha flux. 
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Figure 1. Ha flux per Mq for our synthetic clusters. The shaded contours represent the density of individual realisations of different 
clusters for two methods of filling the IMF: PSS for upper panels and CMDMSM for the lower panels. The lines show the mean fluxes 
from these two methods, the solid black line for PSS and dashed black line f or CMDMSM. The lin es are shown in both single star panels 
and both binary star panels. The black boxes represent the observations of lCalzetti et al ] l|20ld l. the lower points include clusters not 
detected in Ha and the high flux va lue does not in c lude t hem. The black open box re presents the range of possible values derived from 
the Velorum clust er estimated from [jefferies et al ] l|200gf ) and |Pe Marco et al] l|200d l. The diamonds and crosses are based on data of 
lLamb et al. The diamonds represent clusters with upper limits to their mass and the crosses with measured cluster masses. The 

grey triangle indicates the position of the massive HII region NGC604 with values from Eldridgc & Rclafio (2011). The left panels are 
for single-star populations and and the right panels are for binary populations. The linear features in the contours in the panels are due 
to the limited resolution of the stellar models initial masses and time bins used in the synthesis. 



However the observed points at 1O 3 M0 are less conclusive. 
The point without the non-detections lies on the the PSS 
line, while the point including the non-detections lies in be- 
tween the PSS and CMDMSM lines. Thus PSS gives a bet- 
ter fit but a refined CMDMSM scheme that allows a higher 
maximum ma s s for a certain cluster mass may match the 
ICalzetti et all l|20ld ) data. 

An alternative method to discriminate between PSS and 
CMDMSM is to search for individual massive stars that are 
in low-mass clusters. One example is the Wolf-Rayet star 
7- Velorum, the nearest Wolf-Rayet star to the Sun in the 
Galaxy. It is a binary system containing stars that were ini- 



tially 35 and 30 Mp) in a cluster with a total mass of between 
250 and 350M Q (|De Marco et al.ll2000l ; IJefferies et afll2009l : 
lEldridgei r2009'). We have indicated the location of this clus- 
ter in Figure [T] We see that the PSS clusters overlap with 
the parameters of this cluster. The CMDMSM models for 
a single star population do not reach this region. The bi- 
nary CMDMSM models do reach the parameter space for 
7- Velorum. The small number of such models indicates such 
clusters would be rare. This suggests that PSS is more likely 
to be in action in nature although a more relaxed form of 
CMDMSM would also fit the observed data. 

Other more extreme examples of low-mass clusters with 
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a single massive star were observed bv lLamb et al.l |2010t ). 
They observed apparently isolated O stars and found low 
mass clusters associated with these stars. Using the stellar 
and cluster masses derived by lLamb et all J2010) and esti- 
mating the ionising flux for the massive star we have plotted 
their clusters in Figure [T] They are only reproduced by our 
PSS method. 



Testi, 



This agrees with previous stu dies by Testi et al. I dl997T) 
Palla fe Nattal l| 19981 , 1 19991 ) and 'lParker fc Goodwin 



2007) who use similar arguments. iMaschberger fc Clarke 



2008) also made a detailed stu dy of all available information 



and a lso favour PSS. However I Weidner. Kroupa fc Bonnelll 
(2010) performed a similar analysis and found that for low 
mass clusters, below 1OOA/0 PSS is favoured but more mas- 
sive clusters appears to have a CM DMSM relation. The ob- 
servations of lCalzetti et al.l ( |201Q ) do not currently favour 
either PSS or CMDMSM. Here we can only agree that PSS 
occurs in low-mass, up to lOOM© , clusters. For more massive 
clusters, of around lOOOM©, it is difficult to differentiate be- 
tween PSS and CMDMSM. Finally for cluster masses more 
than about 1O 4 M0 the differences are less important. 

Finally we note that our results are in line with those of 
IVillaverde. Cervifio fc Luridianal l|2010bl ) . They suggest that 
for cluster masses below 1O 4 M0 there is a highly asymmet- 
ric scatter of the ionising flux around the mean integrated 
values from standard synthesis models because the single 
most massive star dominates the ionising flux of the cluster. 
This manifests itself in our results by the increased spread 
in Ha flux per Mq at low cluster masses. We note that they 
suggest that PSS is more favoured than CMDMSM. 



3.2 The Ha and FUV in 11HUGS galaxies 

iLee et all (|2009l ). iMeurer et ail (|2009l ) and iBoselli et af] 
1 20091 ) have attempted to gain insight into the IMF by look- 
ing at emission from entire galaxies. They brought together 
Ha observations with far UV cont inuum ob s ervati ons. Here 
we concentrate on the results of ILee et al.l (|2009l ) because 
their set of galaxies are a volume limited sample of 315 
within UMpc. The emission of these two spectral star- 
formation rate indicators are determined by the number of 
stars with masses greater than 20 and 3Mq respectively. 
Therefore measuring the ratio of the two fluxes, or the rela- 
tive star-formation rates measured for the galaxies, gives an 
indication of the number of stars in different mass regimes. 
ILee et all <|2009h found that as the Ha flux decreases the 
Ha/U V ratio decreases so there is more UV flux tha n ex- 
pected. Pfl amm- Altenburg. Weidner fc Kroupal (|2009l ) have 
suggested that this turn down is evidence for IGIMF deter- 
mining the galaxy-wide IMF of these galaxies. Their study 
was based on single star models alone. Here we repeat their 
analysis with binary as well as single star models and also 
our stochastic approach to the star-formation history, with 
stellar clusters forming independently from one another. 

We plot our synthetic galaxies in Figure[2] We see that, 
at star-formation rates above lO _2 M0yr _1 , the spread of 
models is similar but CMDMSM gives a slightly greater 
scatter towards lower values of the Ha/UV ratio. This 
can be more easily seen in Figure [3] where we bin the 
synthetic and observed galaxies with star-formation rates 
above W^M^- 1 by their Ha/UV ratio. CMDMSM has 
a greater range of ratios because of the relative lack of mas- 



sive stars in the total stellar population. We see CMDMSM 
reproduces the lowest ratios at the highest star-formation 
rates. PSS produces much higher ratios. However in this 
model we have assumed no leakage of any ionising photons. 
This would reduce the con tribution from the Ha flux by up 
to around 50 per cent (see lZurita et al.ll2002l . for example). 
This would lead to lower ratio values at high star-formation 
rates for both PSS and CMDMSM. 

To account for the leakage or loss of ionising photons 
from a galaxy or absorption by dust grains we have made 
a simple adjustment to our models. In Figure [3] we have 
modified our synthetic ratio distributions by assuming that 
galaxies lose between and 50 per cent of their ionising pho- 
tons. We take our synthetic populations and smear them 
by this range of possible leakage fraction so that the mean 
leakage is 25 per cent. Even for this modest loss of ionis- 
ing photons the ratio distribution changes the PSS model 
to match the range of observed galaxies. At the same time 
the CMDMSM method has a slightly worse agreement. To 
test the significance of these differences we have used a x 2 
test to compare the observed distributions to the synthetic 
populations. We find that without leakage only CMDMSM 
with single stars is a probable match. However with leakage 
only the CMDMSM single star synthetic population is ruled 
out. 

iFumagalli. da Silva fc KrumholJ (|201ll ) used a leakage 
fraction of 5 per cent and stated that their results were not 
dependent on the amount of leakage. This is because they 
compared the amount of Ha to mean values of Ha flux which 
are less sensitive to leakage than the Ha/UV flux ratio (see 
their figure 2). They found that for leakage fractions up to 
40 per cent their results were unaffected. This is within the 
mean leakage of 25 per cent that we apply to our models. 

For the single star population the greatest difference be- 
tween PSS and CMDMSM is seen in the different paths of 
the mean ratio values versus the star-formation rate deter- 
mined from the UV flux. CMDMSM decreases much sooner 
than PSS at around O.lMoyr^ 1 . However there is a large 
possible range around these mean values in both cases and 
the difference is only approximately la. When we consider 
the binary population we see that the difference between the 
two IMF filling methods is substantially reduced. This is be- 
cause, while the IMF initially leads to fewer massive stars 
in the CMDMSM case, binary interactions, such as merg- 
ing and mass-transfer, increase the number of massive stars 
relative to a single-star population. If we are to distinguish 
between PSS and CMDMSM by means of the downturn in 
this ratio we must repeat the analysis that led to Figure [3] 
for lower star-formation rates. We show the result for star- 
formation rates between 10 -2 and 10 -4 MQ/yr _1 in Figure 
[4] By eye PSS provides a better fit to the observed popu- 
lation than CMDMSM. This is because CMDMSM has an 
extended tail of galaxies towards lower ratios. PSS does not 
have this tail. However including binary stars in the syn- 
thetic galaxies reduces it further in both CMDMSM and 
PSS. Also the tail might not be present in the observed data 
owing to selection biases. A \ 2 test reveals that both PSS 
populations are a probable match to the observed distribu- 
tion. The single-star CMDMSM distribution does not match 
the observed distribution. However our binary CMDMSM 
population produces an equally likely fit to the observed 
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log (SFR(UV) / M Q yr" 1 ) log (SFR(UV) / M Q yr" 1 ) 



Figur e 2. The ratio of SFR measured by Ha and UV fluxes verses the SFR from UV flux. The asterisks are the observations of lLee et al.l 
(2009) while the shaded region show the density of our individual realisations of synthetic galaxies. The thick solid lines indicate the 
mean ratios for the synthetic galaxies and their lcr limits. The dashed lines show the mean ratios for the other IMF filling method with 
the same stellar population. The upper panels are for PSS and the lower panels are for CMDMSM. While the left panels are for a single 
star population and the right panels are for binary populations. Here we assume that a star formation rate of 1M© yr -1 is equivalent 
to a log 10 (F(Ha)/ergs s _1 ) = 41.1 and a UV flux density of log 10 (F(1500A)/ergs s^ 1 A) = 27.85. Linear features are due to limited 
resolution in initial mass, separation and mass ratio parameter space of our binary models. 



data. Our results also show that some ionising photon leak- 
age is required if our PSS models are to match observations. 

iLee et all l|2009t ) noted that their results indicate a 
downturn in the Ha to UV ratio at low star-formation rates. 
Th is could be explained by the IGIMF model pu t forward 
by IPflamm-Altenburg. Weidner fc Kroupal l|2009fl . At first 
comparison of the synthetic and observed galaxies in Fig- 
ures [2] and [3] tempts us to agree with this deduction. This is 
mainly because the spread of the observed galaxies at higher 
star-formation rates is better reproduced by the CMDMSM, 
single star models. The most significant difference between 
PSS and CMDMSM is in the region where the star forma- 
tion rates drop below lO -2 M0yr -1 . All our models are able 
to reproduce the observed galaxies with the lowest ratios 
at low star-formation rates. Therefore it is not possible to 



differentiate between PSS and CMDMSM from these obser- 
vations. Furthermore the inclusion of binaries in stellar pop- 
ulation models means that any difference between PSS and 
CMDMSM is only apparent a t star-formation rates below 
those in the observed sample of I Lee et al.l |2009l 'l. Therefore, 
from the observed distribution of Ha to FUV ratio, it is not 
possible to discriminate between PSS and CMDMSM owing 
to the uncertainties in the importance of binary evolution 
and ionising photon leakage. 

Our conclusions are broadly in line with those of 
iFumagalli. da Silva fc Krumholj ifeoill ). However they com- 
pared PSS models to IGIMF models. The IGIMF models re- 
strict the number of massive stars in the synthetic galaxies 
further because they impose a maximum cluster mass that 
depends on the total star-formation rate. We have only im- 
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Figure 3. The distribution of Ha to UV ratio for observed and synthetic galaxies with star-formation rates between 10 2 and 1 Mqy 1 . 
The red line represents the observed sample of lLee et al . (2009) while the solid line represents the relevant synthetic galaxies from Figure 
[2] The dashed line represents the synthetic observations smeared by a flat leakage of ionising photons distributed between a leakage 
fraction of and 50 per cent. The left panels are for PSS and the right panels are for CMDMSM. The first and third panels are for a 
single star population and the second and fourth panels are for binary populations. 



posed a cluster mass that dependent maximum stellar mass 
and have shown that a CMDMSM alone cannot be ruled 
out. 

An important conclusion to draw from our models (and 
those of iFumagalli. da Silva fc Krumholj|201ll ; IWeisz et al.l 
l2012t ) is that the scatter and variation of the Ha /UV ratio is 
not due to the IMF filling method but it depends more on the 
star-formation history of each individual galaxy. A general 
trend we find is that those systems with less star-formation 
in the last 10 Myr have lower ratios, while those with most of 
the star formation in the last 10 Myr have higher ratios even 
at low mean star-formation rates. This is because the stars 
responsible for Ha emission typically have ages of 10 Myr 
or less. This indicates that any simulation that predicts the 
properties of a sample of galaxies must take into account the 
stochastic nature of star-formation and recognise not only 
that each cluster has its own stellar content but also that 



each cluster has its own age independent of the other clus- 
ters. If there are enough clusters in a galaxy this leads to 
an average stellar population. However if there are only a 
few clusters the appearance of the galaxy-wide stellar pop- 
ulation can be very different from what might be expected 
for a simple stellar population with a smooth star-formation 
history. 

3.3 The importance of binaries 

From our results it is possible to qualitatively demonstrate 
the need to use binary star models. For individual clusters 
binaries seem to have little affect. This is because of the short 
period of 8 Myr we have used to match the observed clusters 
in this case. In our synthetic galaxies, with 100 Myr of star- 
formation, we see that the scatter of the synthetic galaxies is 
reduced slightly if binaries are included and the mean SFR 
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Figure 4. The distribution of Ha to UV ratio for obse rved and synthet ic galaxies with star-formation rates between 10 4 and 10 2 
M@yr~ 1 . The red line represents the observed sample of I Lee et al . (2009) while the solid line represents the relevant synthetic galaxies 
from Figure [2] The dashed line represents the synthetic observations smeared by a flat leakage of ionising photons distribution between 
a leakage fraction of and 50 per cent. The left panels are for PSS and the right panels are for CMDMSM. While the first and third 
panels are for a single star population and the second and fourth panels are for binary populations. 



ratio starts to decrease at lower SFRs. This is more clearly 
shown in Figure [S] in which we compare populations with 
different IMFs and single-star to binary star ratios. We see 
that for a single-star population the ratio begins to drop be- 
tween 10 -2 and 10 -3 Moyr _1 while for binary populations 
this drop begins between 10~ 3 and 1O _4 M0 yr _1 . The bi- 
nary effect makes it difficult to distinguish between the PSS 
and CMDMSM IMF filling methods at any star-formation 
rate. 

Binary evolution affects the observed SFRs because 
through mass transfer between and merging of stars it in- 
creases the number of massive stars at the expense of lower- 
mass stars. We demonstrate this in Figure [5] We see here 
that binary populations typically produce similar Ha/UV 
flux ratios to a single star population when the IMF slope 
is shallower. That is until low star-formation rates at which 
point, for a single cluster with a significant population of bi- 



nary stars we can also expect the apparent IMF to be flatter. 
Furthermore the most massive star in a cluster might not 
have been the most massive star when it formed. Therefore 
interacting binary stars have a strong effect and must be in- 
cluded when attempts are made to determine the IMF from 
observations of stellar systems. 



4 CONCLUSIONS 

We have investigated two uncertainties in population syn- 
thesis. These are how the IMF is filled and the effects 
of interacting binary star evolution. The Ha flux per Mq 
observed in samples of clusters is consistent with PSS of 
the SIMF for clusters around IOOMq because individual 
low-mass clusters with one or two massive OB and WR 
stars, such as the Velorum cluster or those presented by 
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Figure 5. Similar to Figure[2]but here showing the mean ratios 
for single star (solid lines) and binary star (dashed lines) popula- 
tions calculated for PSS with different slopes for the SIMF. 



iLamb et all i|2010r ). provide a strict test to distinguish be- 
tween PSS and CMDMSM. For more massive c luster s com- 
parison with the observations of ICalzetti et al l (l20ld ) is not 
significant enough to rule out CMDMSM. At these masses, 
around 1O 3 M0 and above, it also becomes more difficult to 
to differentiate between PSS and CMDMSM because of the 
blurring effect of binary stars and in addition to the lack of 
conclusive data in this mass range. 

We have also considered the ratio of the Ha to UV 
fluxes in galaxies. Observationally there is a significant scat- 
ter that can be explained by the stochastic nature of the 
star-formation history. We find it difficult to differentiate 
between PSS and CMDMSM. This is because we find some 
evidence that the leakage or loss of ionising photons must 
be considered. In addition, including binary star populations 
makes it difficult to distinguish between the methods for fill- 
ing the IMF. Only single-star CMDMSM populations can be 
ruled out with the observations of galaxies with SFRs below 
lO^Meyr- 1 . 

The ratio of Ha to UV flux for stellar populations, in- 
cluding binary stars, varies less than that for populations 
of single stars. Binaries can merge and mass-transfer can 
produce more massive stars than were present in the initial 
population. Therefore the expected star-formation rates for 
galaxies in which it will be possible to detect differences be- 
tween PSS and CMDMSM are much lower than currently 
observed. Furthermore, because the leakage or loss of ion- 
ising photons from young stellar populations must be con- 
sidered, it becomes even more difficult to discern the IMF- 
filling method from observations of galaxies with low Ha to 
UV ratios. We suggest that it may be more fruitful to find 
galaxies with low overall star-formation rates but with high 
Ha to UV ratios. That is galaxi es th a t are rich in clusters 
similar to those found bv lLamb et al.l (2010). 
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